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We have constructed the quasi-exactly-solvable two-mode bosonic realizations of su{2) and su{l, 1) 
algebra. We derive the relations leading to the conditions for quasi-exact-solvabilty of two-boson 
Hamiltonians by determining a general procedure which maps the Schwinger representations of the 
, sii(2) and su{l, 1) algebras to the Gelfand-Dyson representations respectively. This mapping allows 

' us to study nonlinear quantum optical systems in the framework of quasi-exact-solvabilty. Our 

I approach also leads to a simple construction of special functions of two variables which are the most 

, ^ , appropriate functions to study quasi-probabilities in quantum optics. 

PACS numbers: 03.65.Fd, 03.65.Ca,02.30.Sv 
CN ; INTRODUCTION 

" ^ , The exact solutions of the multi-bosons systems play an important role in quantum mechanics. The Hamiltonians 
■ are employed in quantum optics most often involving multi-bosons. During last decade a great attention has been paid 
t — , to examine different quantum optical models with Hamiltonians given by multi-bosons[l-5]. The quantum optical 
• systems whose Hamiltonians expressed by Casimir operator of the symmetry group can exactly be solved[6]. The 
quantum optical systems with Hamiltonians including nonlinear functions in Lie algebra generators can be analyzed 
by using mainly numerical methods, because Lie algebraic techniques are non-efficient for such systems and most of 
. the other analytical techniques require in general tedious calculations. 

In quantum mechanics there exist physical systems whose finite part of spectrum can be exactly obtained in closed 
forms. They are known as quasi-exactly-solvable(QES) systems[7-ll, 14]. Dating back over fifteen years there has 
I ' been a great deal of interest in QES systems. The one dimensional QES systems based on s/2-a'lgebra have been 
[ classified by Turbiner[7]. The usual approach to the analysis of the one dimensional QES systems is to express 
I the Hamiltonian of the corresponding physical system as linear and bilinear combination of the generators of the SI2- 
algebra. The necessary conditions for the normalizabilty of the wavefunction of the QES systems based on SI2- algebra 
were completely determined in[8]. The QES models, either in the forms of differential equations or of the single boson 
Hamiltonian, have been extended by Dolya and Zaslavskii[12, 13]. 

The aim of this paper is to extend the QES systems to the multi-boson systems and to obtain the solution of 



o 
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^ \ the multi-boson Hamiltonians, in particular, Hamiltonians of the nonlinear quantum optical systems. The formalism 
- - ■ presented here also leads to the construction of special functions with two variables which are appropriate to study 
many quantum optical problems and quasi-probabilty[15, 16]. We think that the QES bosonic systems deserve special 
treatments for they have many applications in physics. The suggested approach can be extended to the matrix 
multi-bosonic systems. 

In the Bosonic Hamiltonian, it is more convenient to use the bosonic representations of the su(2) and su(l,l) 
generators and determine the conditions of quasi-exact solvability. The single boson realization of the su{l, 1) algebra 
has been studied in[12, 13]. In this work we follow a different strategy to obtain the conditions for QES of the 
two-bosonic systems. 

The algebras of the groups SU{2) and S'C/(1,1) have useful appHcations in quantum physics[17-19]. They are 
used to generate the energy spectra while the representation matrices of the group could be used to calculate time 
dependent excitations of the bound states and the scattering states respectively [6]. As usual the group elements and 
their associated algebras can be expressed in terms of bosonic operators. 

This paper is organized as follows. The constructions of the two-mode bosonic realization of the su{2) and su(l, 1) 
algebra are briefly reviewed in section 2. In section 3, we propose the transformation of the boson creation and 
annihilation operators leads to the two class su{2) and two class su{\, 1) one dimensional realizations that is useful to 
study QES systems. The one dimensional realizations of the su{2) and su(l, 1) algebras are presented by similarity 
transformation section 4. We also discuss the some applications of our approach. We demonstrate the procedure 
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presented here can be used to solve a rich family of physical systems. In particular, the Karassiov-Klimov Hamiltonian, 
Hamiltonians of the second and third generation harmonic oscillator systems, Hamiltonian of the quantal systems 
under thermal effect have been solved in section 5. We also demonstrate the relations between QES one-dimensional 
Schrodinger equation and two-mode bosonic Hamiltonians leads to the construction of the two dimensional differential 
equations of the special functions, in this section. The paper ends with a brief discussion and conclusion. 

CONSTRUCTION OF TWO-BOSON HAMILTONIANS: 6u(2) AND su{l, 1) ALGEBRAS 

A convenient way[6] to construct a spectrum generating algebra for systems with a finite number of bound states 
is by introducing a set of boson creation and annihilation operators. We introduce two boson operator, ai and a2, 
obey the usual commutation relations 

[ai,ai] = [ai,aj] = [a2,aj] = 0, [ai,a5^] = [02, aj] = 1. (1) 

The bilinear combinations o^ai, 0^02, ajai and aJo2 generates the group U{2) and a^ai, 0102, aia2 and 03^02 
generates the group SU{1, 1). Let us start by introducing three generators of SU{2), 

J+ = ai"a2, J- = 02 ai, "^0 = ^ (a^ai - 012^02) • (2) 
These are the Schwingcr representation of SU (2) algebra and they satisfy the commutation relations 

[J+,J_] = 2Jo [Jo,J±] = ±J± (3) 
The fourth generator is the total boson number operator 

N = ai + aja2 (4) 
which commutes with the SU (2) generators. The Casimir operator of this structure is given by 

J = J-J+ + Jo{Jo + l) = ^N{N + 2). (5) 
If we denote the eigenvalues of the operator J by 

J = j{j + l) (6) 
the irreducible representations of SU (2) are characterized by the total boson number 

j = N/2 (7) 

where N = 0,1,2,... . 

The Schwinger representation of su{l, 1) algebra can be constructed by considering bosonic realizations of the 
generators: 

K+ = a'^a2, = a2ai, ii'o = ^ (a]'"ai + aja2 + l) . (8) 

They satisfy the commutation relations 



[K+, K_] = -2Ko [Ko, K±] = ±K±. 



(9) 
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The total boson number operator of this algebra is given by 

L = afai — a^a2-, (10) 

and it commutes with the generators of the su{l, 1) algebra. The Casimir operator of this structure can be expressed 
in terms of number operator, such that 

C = -J_J+ + MJo + 1) = i(L + 1){L - 1). (11) 
It is obvious that number operator and Casimir operator commute. If the eigenvalues of the operator C is denoted by 

c = e{e + i) (12) 

the irreducible representations of su{l,l) algebra are characterized by a total boson number 

L=-{2£+l). (13) 

For the exactly solvable case the Hamiltonian have been written in terms of Casimir invariants of the algebra. In 
this case the eigenvahies and eigenfunction of the Hamiltonian can be obtained in the closed form. Our task is now to 
obtain the conditions for QES of the su{2) and su{l, 1) algebra given in (2) and (8). The conditions can be obtained 
by connecting two-mode bosonic realizations of the su{2) and 8^(1, 1) algebra to the sl{2,R) and shiR) algebra, 
respectively. Since Casimir operator J commutes with TV, and C commute with L the linear combinations of the 
generators of each algebra can be diagonalized within the representation [N] and [L] , respectively. The abstract boson 
algebra can be associated with the exactly soluble Schrodinger equations by using the differential operator realizations 
of boson operators. This connection opens the way to an algebraic treatment of a large class of potentials of practical 
interest. 



TRANSFORMATION OF THE BOSONIC OPERATORS 



In the previous section we have summarized the construction of the two-mode bosonic realizations of the su(2) and 
su{l, 1) algebras. In this section we develop a procedure to transform the generators in the Schwinger representation 
to a more suitable representation to determine its relation with the QES systems. This can be done by various 
methods. Here we follow a different strategy to obtain the connection between the two-mode bosonic systems and 
QES systems. Let us introduce the following similarity transformation induced by the operator 



S = (a+)""^»i (14) 

where a is a constant and in order to obtain Gelfand-Dyson representations of the su{2) and su(l, 1) algebra it will 
be constrained to ±1. In general, it will be shown later, it is not necessary to set a = ±1, to construct QES systems. 
This property allow us to study a wide range of physical systems. The operator S acts on the state |ni, 712) as follows. 



n2! 



5|ni,n2) = (02 \ni,n2) = \ — rr |ni,n2 + ani) . (15) 

y (n2 -t- anij! 

Since oi and a2 commute, the transformation of ai and under S can be obtained by writing 02 = e^, with 
[ai,6] = [a+ fo]=0, 



ai(aj)"' 



(16) 
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and transformation of 02 and at as follows 



Sa+S-' = {atr^taia+{a+r^^t^- = at 
Sa2S-'^ = {a+)°"'^''^a2io+)-°"'^''' =a2-aa+ai{a+)-\ (17) 

In a similar manner, from Schwinger representation to the Gelfand-Dyson representation which is suitable to study 
QES systems, the generators of 5^(2) and su{l, 1) algebra in the Schwinger representation can be transformed by 
introducing the following operator: 



The operator T acts on the two-boson state as 



T = ar^"' (18) 



T\ni,n2) = 02"' \ni,n2) = J — — rr |ni,n2 - an\) . (19) 

y (n2 - ani)\ 

Since ai and 02 commute, the transformation of a\ and under S can be obtained by letting 02 = e'^ with 
[ai,c] = [af ,c] = 0, 

TaiT-i = e"'^'*^''iaie-"™^'"^ =ai(a2)-" (20) 
The transformation of 02 and a^ are as follows: 

1 aj J = a2 ajaj = + aaf 0102 . 

J 02-/ = 02 0202 = O2 (21) 

These two transformation leads to the two different su{2) and two different sm(1, 1) realization in one dimensions and 
the corresponding realizations are useful to study QES systems. 

DIFFERENTIAL REALIZATIONS OF THE SU(2) AND SU(1,1) ALGEBRAS 

The realizations (2) and (8) can bo transformed in the form of the one dimensional differential equations in the 
Bargmann-Fock space when the boson operators realized as 

ai = ^, at = X. (22) 

We can obtain two different differential realizations of the su{2) and su(l, 1) algebras, depending on the choice of the 
a, in the equations (14) and (17). 



su{2) Realization 

When the generators (2) of su{2) algebra is transformed by using the similarity transformation operator S in the 
case of a = 1, takes the form: 





= SJ+S-'^ = -{a+ 


fai + N'a+ 


J'_ 


= SJ_S-'^ = oi 




Jo 


— SJqS^^ = o^oi 


N' 


N' 


= SNS-^ = 0+02. 
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These representations are called Gelfand-Dyson representation of the su{2) algebra. Similarly we can easily obtain 
a second realization by using the transformation operator T and a = —1. In this case the realization of su(2) is given 

by 





= TJ+T-'^ 


= at 


J'_ 


= TJ_T-^ 


= -a+{aif + ai{N' 


4 


= TJqT-^ 


+ 

= ajai - — 


N' 


= TNT-'^ : 





(24) 

The difference between the Schwinger and Gelfand-Dyson representation is that while in the first the total number 
of a\ and 02 bosons characterize the the system, in the later it is only the number of 02 bosons that characterize the 
system. According to (8) the representation is characterized by a fixed number 2j. Therefore in the Gelfand-Dyson 
representation, the primed generators can be expressed in terms of one boson operator ai. According to the (7) 
takes the values N' = 2j. The realizations (23) and (24)are the well known generators of the sl{2, R) algebra, in the 
Bargmann-Fock space, which play an important role in the quasi-exact solution of the Schrodinger equation. The 
linear and bilinear combinations of the generators form a second order differential equation and according to the 
Turbiner theorem the linear and bilinear combinations of the generators J- (i = -t-, — , 0) are QES. The basis function 
of the primed generators of the su{2) algebra is the degree of polynomial of order 2j, 

Pn{x) = {x\x\--- ,x^^). (25) 

su(l,l) realizations 



By using the similar arguments given in previous subsection we can obtain two different one dimensional differential 
realization for the su(l,l) algebra in the Bargmann-Fock space. The generators of the su(l,l) algebra under the 
transformation of the T when a = 1 takes the form 





= TK+T-^ 


= {a+fai + iL' + l)a+ 


K'_ 


= TK_T-^ 


= ai 


K 


= TKqT-^ ■ 


+ ^L' + l 
= ajai + 2 


L' 


= TLT-^ = 


aJo2. 



(26) 

The other realization can be obtained by transforming the generators with the transformation operator S and choosing 
a = -1: 



K 


= TK+T-^ 


= at 


K'_ 


= TK^T^^ 


= a+ {aif + {L' + l)ai 


K 




- a^ai + 2 


L' 


= TLT-^ = 


02" 02- 



(27) 

The basis function of these generators are polynomials in x, in the Bargmann-Fock space. In the representations (26) 
and (27) the operator L' characterize the system, L' = —2£ — 1. 

Consequently the (quasi)exact solution of the two-mode bosonic Hamiltonians which include linear and bilinear 
combinations of the su{2) or su(l, 1) algebra can be obtained by a suitable transformation. 

APPLICATIONS 



In this section we discuss the applicability of the method to solve the Hamiltonians of the various physical systems. 
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Karassiov-Klimov Hamiltonian 

The method developed in this article can be used to obtain the solution of the various two-boson Hamiltonians. 
Consider the following family of Karassiov-Klimov [3] Hamiltonians 

H = ujia'^ai + u}2<^2'^2 + Kflj^^ftj + Sa^a^'' (28) 

where < r < s, uii and 102 are frequencies of two independent harmonic oscillator and k and R are coupling constants. 
The Hamiltonian was studied in the context of nonlinear Lie algebras [20], for r = s = 2 and r = 1, s = 2. Let us 
consider the transformation of the Hamiltonian H by the transformation operator S, choosing a = r/s, 

H' = SHS~^ = wiaf ai + LJ2{a2a2 — -a^ai) + Kai^{a2a2 — -a^aiY + Ral (29) 

Using the relation (7) and Bargmann-Fock space realizations (22)of the ai and af the Hamiltonian takes the form 

H' = {coi - -)x^ + 2jU2 + KX%2j - -x-^Y + R^ (30) 
s ax s ax dx^ 

The eigenvalue problem can be written as 

H'Pnix)=EPn{x). (31) 

With the basis function (25) the eigenvalue equation (31) leads to the following recurrence relation 

V T 7)} 

(t^, _ _ _ E + 2ju2)Pn{E) + K{2j - -nyPn+s{E) + R-^^^—^Pn-s{E) = (32) 
and the wave function of the Hamiltonian H can be written as 

V'(x) = 5-i^Pfe(£;)x'= (33) 

fe=0 

The wavcftinction is itself the generating function of the energy polynomials. The eigenvalues are then produced by 
the roots of such polynomials. If the E^ is a root of the polynomial Pk+i{E), the series (22) terminates at > 2j^ 
and Ek belongs to the spectrum of the corresponding Hamiltonian. The eigenvalues are then obtained by finding the 
roots of such polynomials. 

The Hamiltonian (28) have been considered in[20], in the context of the nonlinear algebras, for the specific values 
of s = r = 2 and s = 2,r = 1. The Hamiltonian (23) can be expressed in terms of the generators of su{2) algebra 
when s = r = 2 

H ={loi- lo2)Jo + kJI + RJl + i(u;i + i02)N (34) 
and the transformed Hamiltonian can be cast into a differential operator 

d? d 

H = (R + Kx'^)—^ + x{ux -i02+ 2k(3 - 2j)x^)- 2j(wi + k{1 - 2j)x'^), (35) 

dx'^ dx 

by using the realization given in (23). Our result is coincides with the result given in[20]. The recurrence relations, 
with the basis function (25) is given by 

K{2j - n){2j - n - 1)P„+2(A) + Rn{n - 1)P„-2(A) 

+ {uji{n - 2j) - uj2n - E)Pn{\) = (36) 

The Hamiltonians of the second and the third harmonic generation are the special case of the Karassiov-Klimov 
Hamiltonian and they are given by 

H = Wia^ai + Ci;2aja2 + K(a^^a2 + ajaj) (37a) 
H = uiafai + uj2a^a2 + K{af^a2 + a? 02^) (37b) 

respectively. It is obvious that (28) can be put in the form of the (37a-b) and they can be expressed as one variable 
differential equations as in(30). Their eigenvalues and eigenfunctions can be obtained by using the recurrence relation 
(32) and (32), respectively. 
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Quantal system under thermal effect 

The next example is the Hamiltonian of the quantal system under thermal effect [21]. In order to investigate 
thermal effects in quantum many-particle systems the Hamiltonian can be formulated by using two bosons creation 
and annihilation operators. The Hamiltonians of the symmetric and asymmetric rigid rotators have almost the same 
structure and their algebraic structure is sm(1, 1) algebra. The Hamiltonians of such systems can be expressed as [22]: 

H = fkjj{a'^ai — a^a2) — 17/1(0]'' — 0102) (38) 

Using the realizations (8) of su{l, 1) algebra we can rewrite the Hamiltonian: 

H = hujL- ijh{K+ -K-) (39) 

The Hamiltonian can be expressed as one dimensional differential equation by using the transformation operators S 
and T: 

H' = SHS-^ = huj{-2j - 1) - i^h{{x'^ - 1)^ - 2jx) (40a) 

H' = THT-^ = hu(-2j - 1) - i-yhix - x—^ + 2j—) (40b) 

dx'' ax 

respectively. According to the Turbiner theorem[7] this equation is exactly solvable. 

Two variable difFerential equations of the special functions 

Let us consider the following bosonic operator 

H = (jjiafai + 0)20,2 (12 + aia^aj + 020102 — -^o^ai- (41) 
In terms of the generators of su(l, 1) the Hamiltonian can be expressed as 

^) + ^{oJi-i^2)L + aiJ. 
The transformed Hamiltonian can be expressed as a differential operator 



H = (wi +i^2){Ko - 2) + 2^^^ ~ '^2)^ + + J- - ^Jl (42) 



1 (fi d 
= --^ + (a2 + x{u>i +W2 + "la;))^ - (2fcaia; + 0Ji{2k + 1)), (43) 

which can be written in the form of an eigenvalue equation 

HR{x) = ER{x). (44) 
The solution of (43) can be obtained by introducing the wave function 

i?(a;) = e--f'^(^)'^^t/>(:c) (45) 

where the weight function W{x) is given by 

W{x) = -a2 - (wi + uJ2)x - a-^x^) (46) 

With the wave function given in (45) the Hamiltonian (43) can be transformed in the form of Schrodinger equation 
with the potential 
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V{x) = ^{al - wiiAk + 3) - C02) + {a2{u)i + W2) - ai{2k + l))x + 

i(aia2 + (wi + i02f)x^ + ai{u!i + lo2)x^ + yx^. (47) 

Therefore the bosonic operator given by (41) is related to the anharmonic oscillator potential. It is known that 
the Schrodinger equation with the potential given in(47) is QES. Let us turn our attention to the Hamiltonian (43). 
When we set 

ai =a2 = 0,u}i =W2 = 1/2, E= ^{n-2k-l) (48) 
then the eigenvalue problem can be written as 

1 d n\ ^, . ^ 

-^-r^+x-;---]R(x)=0 (49) 

and the bosonic Hamiltonian takes the form 

H = -{afai + 02 a2 — ala^)- (50) 

It is obvious that differential equation of the two dimensional Hermite polynomials in the Bargmann-Fock space can 
be expressed as: 

+ X,— +X2—]RixuX2)=0. (51) 



2 \ dxidx2 dxi dx 



Our last example is the two-boson realization of the Schrodinger equation with the sextic harmonic oscillator 
potential. Let us consider the following bosonic operator: 

L = {ivi + ai)af ai + {cui — ai)a2 a2 + a+ai a2 + (52) 
(a_ - ^)ai02 - ^{a^aial + a|a2af ) + wi. 

In terms of the generators of the su{l, 1) algebra it can be written as 

H ^ 2LJiJn + a+J+ + a-J- + aiN ~ JqJ-, (53) 
which can be converted to the differential operator 

H = -x-^ + {a-+k + x{2oji + a+x))-^ - 2k{a+ + {ui + a+x)). (54) 
Let us change the variable and redefine the wave function for the eigenvalue problem 

Hij{x) = E^{x),x= (|)^V'(^) = e^'^^'^'^^R{z) (55) 
where the weight function W{z) is given by 

a+ 3 wi l-|-2a_-|-2fc 
^^'^ = l6 " + + 2z • ^^^^ 
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Upon this substitution the potential term in the Hamiltonian can be written as: 




(1 + 2a- + 2fc)(3 + 2a- + 2k) 
4^2 



+ 



2 a+c^ 4 /a+x2 
16 V16/ 



(57) 



This is the radial sextic oscillator potential studied in the literature. 



CONCLUSION 



In this paper we have discussed solution of the two-boson Hamiltonians and applications in physics. It was shown 
that the solution of the two-boson Hamiltonians can be transformed in the form of one- dimensional QES differential 
equation by applying a suitable similarity transformation. This transformation also leads to the connection between 
one and two-dimensional special functions of the physics. 

The method given here is useful to study nonlinear quantum optical systems. The range of the Hamiltonian can 
be extended by the Bugoliov transformation of the boson operators. It is expected that this work will leads to the 
construction of the multi-boson QES Hamiltonians and their extensions to the matrix Hamiltonians. We have been 
presented a first step towards an extension of the quasi-exact solution of the multi-bosonic Hamiltonians. 
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